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Abstract 



u 

^^ ' The redundancy allocation problem is formulated minimizing the design 

cost for a series-parallel system with multiple component choices whereas 
ensuring a given system reliability level. The obtained model is a nonlinear 
integer programming problem with a non linear, non separable constraint. 
We propose an algebraic method, based on Grobner bases, to obtain the ex- 
act solution of the problem. In addition, we provide a closed form for the 
required Grobner bases, avoiding the bottleneck associated with the compu- 
^ . tation, and promising computational results. 

I> 

o. 

^ ■ 1 Introduction 

Slj ' System reliability is considered an important measure in the engineering design 

^ ■ process. A series system is like a chain composed of links, each of them repre- 

senting a subsystem. The failure of one of these components means the failure 
of the whole system. In order to avoid this, it is usual to use redundant compo- 
k>( ' nents in parallel to guarantee a certain level of rehabihty. These systems are called 



^ . series-parallel systems. 



Determining the optimal number of components in each subsystem is the so 
called reliability optimization problem. Two different approaches are usual: 

• maximize system reliability subject to system budget constraint, or 

• minimize system cost subject to a required level of reliability 

Both problems are nonlinear integer programming problems, and they are NP-hard 
|[5l . There are very few papers looking for their exact solutions, due to the difficulty 
of the problems. Those works use essentially Dynamic Programming [14], branch 
and bound methods 19], or Lagrangian relaxation [11], among others techniques. 
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On the contrary, in the Uterature there are many heuristics and metaheuristic algo- 
rithms; such as those based in Genetic Algorithms ||6l, Tabu Search iflOl or Ant 
Colony Optimization (T], among others. 

In this paper we study the exact solution of one of the versions of the prob- 
lem that minimizes the cost function of the chosen design, subject to a non linear 
constraint which describes the reliability of the considered system. For a fixed 
subsystem its inner components can be considered equal, as in [9], or different, as 
in ifTTl . If the components are equal the reliability function is separable and con- 
vex, and the problem can be reduced to a linear knapsack problem [9]. In the case 
of multiple component choices the reliabihty function is no longer separable. In 
|[T4l . the solution is found using dynamic programming methods. That approach 
presents two stages; in the first one the problem is restricted to each subsystem, 
with a level of reliability. Under this assumption the reliability function is separa- 
ble, and the optimization problem can be reduced to a knapsack problem. Then the 
reliability levels of the subsystems are determined by a new dynamic programming 
process. 

The solution method shown in ifTTTl uses an algorithm based on Lagrangian 
relaxations over two linear relaxations of the original problem. The first relaxation 
consists of deleting the non-linear reliability function, and adding certain linear 
constraints, one for each subsystem. The second relaxation assumes that the same 
type of component is going to be used in every subsystem, so that the problem has 
the form as in |[9l. 

Mainly, the algorithm of [H] is a what their authors called a cut and partition 
scheme (a geometric branch and bound). The solution space is partitioned in boxes, 
which are divided and discarded for certain conditions. The cuts are built from the 
best bound feasible solution of some Lagrangian relaxations. Such bounds allow 
to remove certain boxes depending on the improvement with respect to the current 
best point. 

We address the problem via a different approach based on Grobner bases. As 
introduction on this subject, we recommend the text books lUl, ||3l and lHJ. 

Grobner bases were applied to Integer Linear Programming, by the first time, 
in IlZJ. Later, Tayur et al. [13| introduced a new application framework, which 
solves nonlinear integer programming problems, with a linear objective function. 
This is exactly our framework, as in |@]. 

First, we consider a relaxed integer programming problem where all the restric- 
tions are linear. Then we find the solution of the relaxed problem by computing a 
test set. By using the so called reverse test set, we can solve the complete prob- 
lem, generating paths from the solution of the relaxed problem to a solution of the 
complete one. These paths increase the cost function at each step. 

A test set for a linear integer programming problem is a set of directions that 



can be used to design descending algorithms with respect to a linear cost function. 
A test set can be computed from a Grobner basis of the toric ideal associated with 
the linear restrictions, with respect to an order given by the cost function. 

One of the main tasks in the process described before is usually the calculation 
of the Grobner basis. We construct a linear programming problem from the original 
one, removing the reliability function, and adding a new linear restriction. This 
constraint is obtained computing a feasible solution with a greedy algorithm. For 
the relaxed linear programming problem obtained in this way we explicitly give the 
associated Grobner basis, and so the test set to solve the main problem. We point 
out here that any Grobner basis computation is done using closed formulas, thus 
avoiding the hard computation burden of reduction algorithms to compute Grobner 
bases. 

The organization of the paper is as follows. In Section |2] we introduce the no- 
tation and describe the model of a parallel-series system with multiple component 
choices. In Section [3l a greedy algorithm is described to compute a feasible point. 
Section |4] is devoted to a brief introduction to the essential facts about Grobner 
bases. Section \5\ contains the main result about the closed formula for the test 
set of the integer linear problem. Our computational experiments are reported in 
Section |6l Finally we draw some concluding remarks in Section |7J 

2 General model 

In order to formulate the problem, some notation is first introduced. 

• n number of subsystems. 

• kj number of different types of available components for the i-th subsystem, 
/ = \,. . . ,n. 

• rjj reliability of the j-th component for the j-th subsystem, i = l,...,n, j = 

1, . . . ,K,-. 

• Cjj cost of the j-th component for the i-th subsystem, i=l,... ,n, j = l,...,ki. 

• kjjUij lower/upper bounds of number of j components for the i-th subsystem, 
/= l,...,n, j = l,...,ki. 

• Rq admissible level of reliability of the whole system. 

• Xij number of j components used in the i-th subsystem, / = !,...,«, j = 
1, . . . ,K,-. 
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Figure 1 : A series-parallel system with multiple choice components 



In our model, some assumptions are considered: 

• Components have two states: working or failed. 

• The reliability of each component is known and is deterministic. 

• Failure of individual components are independent. 

• Failed components do not damage other components or the system, and they 
are not repaired. 

This model is illustrated in Figure [T] It is a system with n subsystems with the 
notation introduced before. The optimization problem can be formulated as: 

(RP) min ILiIy=iCO-%- 



s. t. 



R{x) > Ro, 




i!j=i^ij>h 


/ = 1,.. . ,n, 


0<lij <Xij <Uij, 


/= \,...,n, 




j= l,...,ki, 


Xij G Z+ 


for all i,j, 



(1) 



where /?(x) =nLi(l ~n/'=i(l ~^ijy'')- The first « linear inequalities in (O assert 
that each subsystem must have, at least, one component. 

As usual, we can make a change of variables ytj = xij — lij, so that we can 
assume Z,y = 0. This does not alter the equations of (RP), and some of the last 
equations can be redundant. Hence it can be assumed /,y = without loss of gen- 



erality, and: 

(RP) min i:U,l!]^,CijXij 
s. t. 

R{x)>Ro, 

l!^'=iXij>l, i=l,...,n. 
< Xij <Uij, i = I,. .. ,n, 

J ^ 1 , . . . , K,- , 

Xij G Z+ for all i,j. 

A feasible solution is sometimes called a reliable solution because it ensures a 
reliability greater than or equal to Rq. 

3 Computing a reliable system with a greedy procedure 

The main step used in the algebraic algorithm described in this article is to consider 
an integer linear programming problem (LRP), relaxed from the original problem 
(RP). There is only one nonlinear constraint in (RP), the equation which ensures 
the reliability of the whole system. Removing the nonlinear constraint we get an 
integer linear programming problem: 



k, 



(LRPl) min rULj=iCijX, 



'J 



s. t. 



-k, 



L/=1-X0>1' i=l,...,n. 
< Xij <Uij, i = I,.. . ,n, 

J ^ l, . . . ,Ki, 

Xij G Z_|_ for all /, j. 

In our solution technique, we start from the solution of the linear programming 
problem {LRPl), and following the directions given by the test set of (LRPl) we 
follow a descent path to the solution of the complete problem {RP). If the linear 
relaxation is too weak, the paths to be followed to get to the optimal solution of 
{RP) will be very long, and the number of points to be processed is huge. 

To avoid this problem we add a new linear equation. We need a feasible point 
3'" of {RP), and there are lots of heuristic methods to obtain such a point. In our 
case, we use a greedy algorithm similar to [9| or [11]. 

At the beginning of the greedy algorithm, 3^" describes the system with the max- 
imum number of components of every type. If the reliability of that system is less 
than Rq, the problem is unfeasible. We consider / the set of all pairs {i,j), which 
describes the j-th component for the i-th subsystem. For each {i,j), we calculate 



the rate tjj = _j^ (i-r-) between cost and reliability, and order / non increasingly 
by these rates (ties are solved by lex order, for example). For the first index (/o, Jo) 
in /, we subtract components of type (/o, jo) from y^ until it is non reliable, there 
is no such component or the io-th subsystem is empty. If the solution obtained by 
this process is non reliable, or the iQ-th subsystem is empty, one (ioJo) component 
is added. Then we take the next index in the set / and repeat the procedure, until 
the index set / has been completely processed. 



Data: r,y, vector c 
Result: )»" feasible point 

y = {un,---,UnkJ 

t = 



_£u 



£nlq. 



log(l-ru)'' ■• ' -log(l-r„t„)^ 
/ = {(l,l),...,(l,^l)...,(«,^„)} 

Order / non increasingly by tij 
forall the (/, j) G / do 

Reliable=TRUE 

S ubsy stemNonEmpty=TRUE 

while Reliable and SubsystemNonEmpty and y*' > do 

ifl-t3'li<lthen 

I SubsystemNonEmpty=FALSE 
end 
if/?(/)</?othen 

I Reliable=FALSE 
end 
if Reliable=FALSE or SubsystemNonEmpty=FALSE then 



^1} 



--y. 



0+1 



end 



end 



end 



Algorithm 1: Greedy algorithm 



Using the above greedy algorithm, we obtain a feasible point j", with a cost 
Hij^ijyh' — ^°- The optimal solution of (RP) has a cost less than or equal to c", so 
we can add to (RP) a valid inequality stating this condition and the problem has an 



equivalent fomi: 










{RP) 


min 

s. t. 


Lr=l L/=l '^iF'j 

R{x)>Ro, 










I'=i^/;>1> 


/ = !,.. 


.,n, 






0<Xij <Uij, 


/ = !,.. 
j = h- 


.,n, 



L;=i L/=i CijXij < c , 

Xij € Z+ for all /,y. 

From this formulation, we have the new integer linear problem 

(LRP) min i:Ul!j=iCijXij 
s. t. 



i:}=i^ij>h 






/ = !,...,«, 


0<Xij <Uij, 








T.i=ilLj=iCijXij 


< 


c°, 




Xij € Z+ 






for all i,j. 



4 A review on integer programming and Grobner bases 

In this section, we recall the concepts and algorithms used to solve Integer Linear 
Programming problems from an algebraic point of view, and the walk back proce- 
dure for nonlinear integer programming problems based on test sets. To this end, 
we have followed 1.12.1 and 1.13.1. 



4.1 Grobner bases 

Denote by /:[jir] =k[xi,... ,xn] the ring of polynomial with coefficients in a field k. 
In our case, k will be M. The ideal generated by a subset J^ C k[x] is the set (^) 
consisting of all linear combinations: 

{.^) = {hifi + --- + hrfr : /l,...,/,G^,/ll,...,/l,€^M}. 
A term order on N^ is a total order ~< satisfying the following properties: 
• -< is compatible with sums, i.e., a ^ j8 =^ a + 7 ^ j8 + 7, for all a, j8 , 7 G N^. 



• ^ is a well-ordering, i.e., -< a for all a G N^, a 7^ 0. 

If we fix a term order -<, then every non zero polynomial / has a unique initial 
term in_<(/) = ojc". It is the monomial ax"" where a is the largest term appearing 
in / for the term order -<. We are particularly interested in two term orders: 

1. The lexicographic order <iex. For every a,/3 € N^, we say a >iex jS if, in 
the vector difference a — j8 € iJ^ , the leftmost nonzero entry is positive. 

2. The vector induced order <c. We consider a vector c G N'^. Given a,j8 € 
W , we say a >c j8 if 

c'a > c'p or c'a = c'j8, and a >iex j8. 

For example, consider the polinomial / = 6x\x\xt, + 7x3 — 5x\ + '^x\x\ and the 
vector c = {3,2,2)' . Then 

in<,ex(/) = -^AM<c{f) = ^AA- 

Of course, we can reorder the variables x,, and get a new term order. In general, 
the notation <c means a term order which respect the partial order defined by the 
vector c and then a tie-break term order, so if we change the lexicographic order in 
the definition on <c, we get another vector induced order. 

Suppose that / is an ideal in k[x], and ^ is a given term order. Then its initial 
ideal is the ideal generated by the initial terms of the polynomials in J: 

in^(7) = (in^(/) : /G7). 

A finite subset ^ of 7 is a Grobner basis with respect to the term order -< if the 
initial terms of the elements in ^ suffice to generate the initial ideal: 

in^(/) = (in^(^) : ge<^). 

Fixed an ideal and a term order, a Grobner basis is not unique. Adding two more 
conditions, the uniqueness is guaranteed. The reduced Grobner basis of / with 
respect to ^ is a Grobner basis ^_^ of J such that: 

• m^{gi) has unit coefficient for each gi G ^_<. 

• For each gi G ^^, no monomial in gi hes in (in_<(^_<\^;)). 
Every ideal J has a unique reduced Grobner basis for each term order. 



4.2 Test set 

Consider a linear programming problem: 

LP{b) min c'x 
s. t. 

Ax = b, 

where A£Z'''^^,b eZ'',c eR^. The notation LP{b) denotes the Unear program- 
ming problem with the right-hand-side restrictions fixed to b. When we write (LP), 
we note all integer programming problems, obtained by varying the right-hand-side 
vector b, fixing A and the cost function c. Consider the map n : N^ — ;■ 7j^ defined 
by n{x) =Ax. Given a vector b € Z'', the set n^^{b) = {« G N^ : n{u) = b} is 
the fiber of (LP) over b. 

We group points in N^ according to increasing cost value dx, and refine this 
order to a total order <c breaking ties among points with the same cost value by 
adopting some term order (lexicographic, for example, as defined in the previous 
section). It is the vector induced order. 

A set G<c C Z^ is a test set for the family of integer problems (LP) with respect 
to the matrix A and the order <c if 

• for each nonoptimal point a in each fiber of (LP), there exists g € G<^ such 
that a — g is a feasible solution in the same fiber and a — g <c cc, 

• for the optimal point j3 in a fiber of (LP), (5 — g is unfeasible for every g G 

A test set for (LP) gives an obvious algorithm to solve an integer program, provided 
we know a feasible solution to this problem. At every step of this algorithm, we 
have two different cases: 

• There exists an element in the test set which, when subtracted from the cur- 
rent point, yields an improved point. We are then in a nonoptimal point, but 
we get a better one. 

• There will not exist such an element in the set, so we are in the optimum of 
the fiber. 

4.3 Toric ideal 

We define Ia the toric ideal associated with A as 

Ia = {x"-xI^ : Aa=Aj3,a,i3GN^). 



Given an integral vector 7 G ll^ , we can write it uniquely as 7 = 7+ — 7 , where 
7+, 7^ € N^ and have disjoint supports. It is well known ( |[T2l ) that 

The relationship between the previous concepts is that the reduced Grobner 
basis ^<^ of /a with respect to the order <c allows us to compute a uniquely defined 
minimal test set G<^ for (LP). The reduced Grobner basis is formed by binomials 

(S<^ = {jc«' - jc'^' , r = 1 , 2, . . . , r} , with in<^ (jc"' - x^' ) = jc"' , 

and then the test set is expressed as 

G<, = {a,-A-,' = l,2,...,r}. 

4.4 Walk back procedure 

Basically the walk back procedure gives an algorithm which computes the optimum 
for a nonlinear integer programming problem under some conditions. The integer 
programming problem {RP) introduced in Section[2]is not linear. It has a nonlinear 
constraint (the reliability condition), while the rest of the restrictions are linear and 
the cost function is also linear. These are the conditions required to use the walk 
back procedure, introduced in ||T31 . In Algorithm |2] it is used the directed graph 
defined by the Grobner basis over the feasible points, but directions are reversed in 
the skeleton. In each step, elements w = a +^ in the reverse skeleton are computed, 
where Aa = and g is an element in the Grobner basis. 

In general. Algorithm |2j uses the following notation. We denote by {RP) the 
entire non linear integer programming problem, (LRP) the relaxed linear integer 
programming problem which arises from [RP). Let j8 be the optimum of (LRP). If 
j8 is feasible for (RP), then it is the solution to (RP). If it is non feasible, then the 
reverse skeleton is needed. 

Let P{cc) denote the path, in the directed graph (reversed) of the linear integer 
programming Problem (LRP), from the optimum j8 for (LRP) to a feasible point a 
for (LP). There is always one. Any solution of (RP) is feasible for (LRP), so the 
objective is to find such a path, in an ordered way. In each reversed step the cost 
function increases, so the minimum cost feasible points for (RP) are found first. 

5 The test set for the relaxed linear problem 

Once the (LRP) problem is reinforced by means of the linear constraint that comes 
after a feasible solution of (RP) is found by the greedy algorithm, the relaxed linear 
problem is: 
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Data: Matrix A, vectors b,c, non-linear restrictions 
Result: Optimum 

^ = {Pili)} 

/ = greedyiRP) 

y = {/} 

§f<^ = groebner(/A) with respect to <c 
j3 = optimum for relaxed LRP 
repeat 

forall the P{a) G ^ do 
forall the g € ^<, do 
w = a + g 

if w is a feasible point of (LRP) then 
if w is feasible for (RP) then 
Y = YU{w} 
Prune P{w) 
ise 

if J <c w/or 5'ome j € F then 
I Prune /'(w) 
end 

^ = ^U{P{w)} 
end 
end 
end 
end 

Delete P{a) from ^ 
until all paths in HP are pruned 
Optimum = Select minimum <c element from Y 



Algorithm 2: Walk back procedure 
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(LRP) min rLil!j=i CijXtj 
s. t. 

0<Xij < Uij, 

Xij G £j^ 



7 = 1 , . . . , Ki , 



for all ij. 



Each inequality must be converted to an equality, so we must introduce a new 
slack variable for each inequality: 

{LRP) min ILi l''=i Qyx^y 
s.t. 

Xij + tij — M,y, 



/= !,...,«, 



LLli:}UiC^;^^7 + ^^ 



X;y G ij-i- 



7 = 1, . . . ,/c,-, 
for all /, j. 



If we put N = ki + . .. + kn and 



^nxW 



0...0 1...1 



K \ 
0...0 



\0...0 0...0 ... 1...1/ 
the restrictions in matrix form can be written as 



D —In OnxW OnxT 

In ^Nxn In 0/vxi 

iCixW Olxn OixW 1 



"nxl 
^/Vxl 

V ^ / 



Ml 



V-7 
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where 

Clx/V 



Cll 



Clh 



Cn\ 



/ x„ \ 




( ''' \ 


Xlki 


,t = 


hk, 


X„i 




tn\ 


\ ^nK ) 




\ tnk„ J 



■■■ Cnk„ ) 




i= \,...,n, 



and 1„ denotes the n- vector with all the componentes equal to 1. We can assume 
that, for each / = l,...,n, the costs c,y are ordered in descending order: ciq > Cip if 
q < p. Let 

Z = {xn, ■ ■ ■ ,Xnk„,di,. . . ,d„,tu,- ■ ■ ,tnk„,b) = {x,d,t,b), 

and consider the following set of binomials in k[z]- 

iox i = \,. .. ,n,k = \,. .. ,ki,\ < q < p <ki. Let > be a term order in k[z] such 
that X > d > t > b. Within each block, the variables are sorted lexicographically 
as follows: 

;Cll > • • ■ > Xu, > X21 > • • • > Xnk„,til > ■■■> tiki >t21> ■■■> t„k„,di > ■■■>dn. 

Theorem 1. The set ^ is the reduced Grobner basis of the toric ideal I a with 
respect to the term orden >. Moreover, ^ is the reduced Grobner basis with respect 
to the order <c induced by the cost vector c. 

Proof. The proof follows similar steps and notation that (IT. Thm. 4]. First of all, 
the set ^ is a subset of Ia, because all the binomials z" — z'^ in ^ verify Aa = Aj8. 
The initial term of every binomial in ^ with respect to > is the underlined term. 
It is enough to show that for every binomial z" —Z^ &Ia, with initial term z", there 
is some g^'S whose initial term divides z". By definition of toric ideal, z"' — '^ € 
Ia if and only if a - j8 G /: = {j € Z' : Ay = 0},^ = « + 2A^ + 1. We denote 
an element y in K by y = {yx,yd,yt,yb) to indicate the correspondence between 
components of j and the columns of A. In addition, we denote the components of 
yx by (Xii , . . . ,Xnk„), and similarly for the others. We classify the elements in K in 
the following manner: 
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1. Let^i = {y £K : yx = 0}. Now j G Ki if and only if {yd,yt,yh) G '^'",si 
n+N+\ belongs to the lattice S' = {w e H^ : A'w = 0} where 



A': 



But 5" = since A' is a non singular matrix. Therefore ^i = 0. This im- 
plies that there are no binomials of the form z" — z^ that do not contain the 
variables x,y. 

2. 'L&iK2 = {y^K : y^ = 0}. Again y eK2 Hand only if {yx,yt,yh) eZ'^,S2 = 
IN + 1 belongs to the lattice 5"' = {w € Z'^ : A"w = 0}, where 



A": 



Let Xiq be the left most nonzero component of y^. We may assume that 
Xiq > since S is the set of integer points in a vector space which implies 
that it contains the negative of every element in it. The /-th row of matrix D 
in A" implies that there exists some p> q such that Xip < 0. Therefore, 

Xiq divides z^ and x,-^ divides z^ ■ 

Consider now the rows given by the block ( /a? I^ ) in A". These rows 
imply that Tiq = —Xiq < and Tip = —Xip > 0. Therefore, 




Xiqtip divides z^ and Xiptiq divides z^ 



The initial term of z^ — z^ with respect to > is z-'' since Xiq divides z^ 
and Xiq is the greatest variable that appears in this binomial. But this implies 
that the initial term of Xiqtip — Xiptiqb'^'i^''''i' € ^ divides the initial term of 
Z^ — Z^ ■ Therefore, the initial term of all binomials associated with K2 is 
divisible by the initial term of an element in ^. 

3. Consider now a general element in 5 = {j € Z* : Ay = 0},s = n + 2N + 1 , 
with no variables restricted to be zero. By the previous cases we may assume 
yx / 0,yd 7^ 0. Let D,- be the first nonzero component of y^- As before, we 
may assume that D, > 0. Therefore, 

di divides z-'' ■ 
14 



Then there exists Xit > 0, so 

Xiicdi divides z^ ■ 

Similarly, there exists Tik < and 

tik divides z^ , 

so the initial term of z^ — 7? is z^ because di > tik. This initial term is 
divisible by Xi^di, which is the initial term of Xikd, — titb'^"'. Therefore 

in<(/^) = (in<(l^)), 

which proves that 1^ is a Grobner basis of Ia with respecto to the term order 
<. Clearly, it is reduced. 

Moreover, with respect to the term order <c, 

in<, (xikdi - tik¥''' ) = Xiudi, 

because the weight of the first monomial is equal to c,i > 0, and the weight of the 
second monomial is equal to zero. Similarly, 

in^^X^iqtip -Xiptigb"-"^''"') =Xiqtip, l<q<p<ki, 

because the weight of the first monomial is Cj^ > cip, which is the weight of the 
second monomial. If Ciq = Cip, the tie is broken with the lexicographical order 

D D 

The above theorem gives a reduced Grobner basis with respect to the term or- 
der induced by the objective function of {LR). On the contrary, [13 Thm. 4] only 
provides a Grobner basis with respect to a lexicographical order, and not with re- 
spect to the term order needed for the computation of the test set. Therefore, in 
order for that Grobner basis to be applied to solve their problem one more compu- 
tational step is required whereas our construction gives directly the answer with its 
consequent saving. 

6 Computational results 

The previous construction of the test set is used in our computational experi- 
ments. In order to gain some insights of its efficiency, if a program like 4ti2 
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(| 4ti2 team(2008)[ ) were used to do the computation of the test set, a simple con- 



figuration of 10 subsystems with 3 components would take more than 60 minutes. 
Therefore, it is very important to apply the result in Theorem [T] to be able to con- 
struct the test set. 

Our algorithm has been coded in M ATLAB and run on a AMD Opteron 252 (2.6 
GHz) with 5 GB RAM. For Table [H all the data in the test problems are randomly 
generated from uniform distributions, with Zj-y = 0, M,y = 4 and r,y € [0.99,0.998], 
as in [llj. The linear cost function Y!i=i L/'=i (^ij^'j ^^^ values ctj G [10,20]. 

The number n is the number of subsystems, and k is the number of different 
components in each subsystem. 

Table 1: /?o = 0.90,r,-y G [0.99,0.998] 



n 


k 


Nodes 


Iter. R-S 


Avg. CPU time (s) 


10 


2 





696 


0.0 


10 


3 





5797 


0.0 


10 


5 





26427 


0.0 


15 


2 





15184 


0.0 


15 


3 


0.4 


85103 


0.1 


20 


2 


7041 


294747 


276.0 



The average CPU time, and the average number of generated nodes during the 
algorithm has been obtained by running the program for 10 instances. 

The column "Iter. R-S" contains the number of iterations according to the re- 
sults given in fTT, Table 1]. Comparing our results with [11], not only the CPU 
time is improved, but also the effort measured by the number of processed nodes 
by the walk back procedure is less than the number of iterations in [11]. We also 
point out that the iterations in lITTI compute two Lagrangian discrete relaxations 
and their corresponding solutions for the best value, each time. After that, to dis- 
card remaining boxes in their branch and bound tree, reliability of that solution 
is needed. In our method in each iteration we only compute a node by adding a 
vector, and then compute its reliability. 

In order to better illustrate the results, new tests have been done with the ad- 
ditional hypothesis that a greater reliability in a component implies a greater cost. 
Note that if there is no correlation between cost and reliability of a component (as 
in nil), then it is likely that certain components are not going to be used, because 
only more reliable components are going to be chosen regardless of their cost. 
Hence the dimensionaUty of the problem is artificially reduced. The results of this 
more realistic case appear in Table |2] It is clear the increasing computational effort 
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Table 2: /?o = 0.90, r/y e [0.99,0.998] (ordered) 



n 


k 


Nodes 


Avg. CPU time (s) 


10 


2 





0.0 


10 


3 





0.0 


10 


5 





0.0 


15 


2 


45 


0.2 


15 


3 


661 


4.4 


15 


4 


28023 


10571 


17 


2 


12578 


1355 



showed by rows n = 15,^ = 3 and n= 15,^ = 4. From the above, we conclude that 
the computational experiments for this model should be done with this additional 
hypothesis of correlation between cost and reliability of each component. 

We also note that the algorithm is very sensitive to changes in the value of 
the reliability parameters r,y. For example, for less reliable components, rij € 
[0.980,0.990] and the same value, Rq = 0.90, for the overall reliability, we have 
obtained the results in Table |3] The reader may observe that the system sizes that 

Table 3: Rq = 0.90, r;y G [0.98,0.99] (ordered) 



n 


k 


Nodes 


Avg. CPU time (s) 


6 


4 


14 





6 


5 


39 


0.1 


7 


4 


1186 


7.4 


7 


5 


5662 


140 


8 


4 


46709 


7010 



can be solved are smaller. However, we have to point out that an exact solution has 
been found in all the examples. An interesting remark is that the elapsed time is 
significantly reduced if the algorithm is stopped with the first best point found in 
the walk back procedure. Obviously, this approach does not guarantee optimality 
but it gives very accurate approximations. From this observation, we think that a 
promising open field is the combination of this technique with heuristic methods 
to get a good approximation of the optimal solution. 
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7 Conclusion 

We have presented in this paper an exact method for solving a nonhnear integer 
programming problem arising from the design of series-parallel reliability systems. 
The method is based on the construction of a test set of an integer linear problem 
through the theory of Grobner bases. We provide an explicit formula of the test 
set, avoiding the high cost of this computation. Computational tests show that this 
approach improves existing methods already applied for this problem. 

This paper deepens the challenge given in IT31 to yield efficient algorithms for 
problems in integer problems based on attractive bases. 
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